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The light-harvesting is a problem of long interest. It becomes active again in recent years stim- 
ulated by suggestions of quantum effects in energy transport. Recent experiments found evidence 
that BChIa 1 and BChla 6 are the first to be excited in the Fenna-Matthews-Olson(FMO) protein, 
theoretical studies, however, are mostly restricted to consider the exciton in BChla 1 initially. In 
this paper, we study the energy transport in the FMO complex by taking different initial states into 
account. Optimizations are performed for the decoherence rates as to maximal transport efficiency. 
Dependence of the energy transfer efficiency on the initial states is given and discussed. Effects of 
fluctuations in the site energies and couplings are also examined. 
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I. INTRODUCTION 

The knowledge of photosynthetic processes in bacte- 
ria and plants may provide us deep insights into design- 
ing efhciently and robustly artificial light harvesting sys- 
tems, where the initial stages in the conversion of solar 
energy into chemical and other forms of energy can be 
described by exciton dynamics with trapping and deco- 
herence. Recent experimental and computational stud- 
ies suggest that environmental noise can assist the exci- 
ton transport and can be optimized for maximal energy 
transfer efficiency |ll-[l3| . In particular, special attention 
has been focused on the so-called Fenna-Matthews-Olson 
(FMO) complex that promotes energy transfer from the 
main light-harvesting complex towards the chemical re- 
action center in green sulfur bacteria. 

As a small protein in green sulfur bacterium, the FMO 
complex is a highly efficient energy transfer wire that 
connects chlorosomes(i.e., light collecting pigment ar- 
rays) to photosynthetic reaction centers. In structure, 
the FMO protein is a trimer, each monomer contains 
seven bacteriochlorophyll-a (BChla) molecules embed- 
ded within a protein scaffold and is believed to function 
independently [Ij]. We shall refer to the individual BChla 
molecules as sites in this paper. Recent studies have de- 
termined the orientation of the FMO complex within the 
inter-membrane re gion between the chlorosome antenna 
and reaction center [lj,[l3, evidence that photosynthetic 
systems exhibit quantum coherence at ambient tempera- 
ture has been found, but the significance of these discov- 
eries for biological energy transfer remains unclear. It is 
quite certain that the reaction center is strongly coupled 
to BChla 3 and that the excitation energy enters an FMO 
monomer from the chlorosomes via BChla 1 or BChla 6. 

Although a preliminary understanding of the complex 
dynamics of excitation transfer has been made, the the- 
ory is still in its infancy. For example, most studies are 
restricted to one exciton limit and the initial exciton is 
supposed in only one of the seven sites. This paper ad- 



dresses two open questions concerning the relevance of 
excitonic coherence to photosynthetic energy transfer. 
First, we present a theoretical framework which eluci- 
dates how decoherence could assist the excitation trans- 
fer, we optimize the decoherence rates for maximal trans- 
fer efficiency and study the robustness of excitation trans- 
fer against the fluctuation of the on-site energy. Second, 
we present a study on the coherence in the initially ex- 
cited states by showing the dependence of the transfer 
efficiency on the initial states. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model to describe the FMO complex and 
the Markovian master equation for the dynamics of the 
exciton transport. In Sec. Ill, we optimize the deco- 
herence rates for maximal excitation energy transfer ef- 
ficiency and explore how the efficiency depends on the 
initial state of the FMO. The fluctuations in the local 
on-site energy and in the site-site couplings affect the 
transfer efficiency, this effect is studied in Sec. IV. Fi- 
nally, we conclude our results in Sec. V. 



II. MODEL 

In the following, we will briefly introduce the model 
Hamiltonian that we use to describe the FMO complex. 
For each site we take two electronic states into account. 
Since we restrict ourself to a single electronic excitation, 
the so-called Frenkel exciton Hamiltonian, 



H = J2Ej\j){j\+ Y. J^A\^){j\ + h.c), (1) 

j=l i>3=l 

can describe the reversible dynamics of the electronic de- 
grees of freedom. Here \j) represents the state where only 
the j-th site is excited and all other sites are in their 
electronic ground state. Ej is the on-site energy of site 
j, and Jij denotes the excitonic couplin g be tween sites i 



and j. In the site basis, we follow Ref.|15j and employ 



the Hamiltonian matrix elements (in units of cm ^) in 
the remainder of this paper, 
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Here the zero energy has been shifted by 12230 cm^^ 
for all sites, corresponding to a wavelength of ^ SOOnm. 
We note that in units of ?i = 1, we have 1 ps^^=5.3 
cm~^. Then by dividing Jij and Ej by 5.3, all elements 
of the Hamiltonian are rescaled in units of ps^^. We can 
find from the Hamiltonian H that in the Fenna-Matthew- 
Olson complex (FMO), there are two dominating exciton 
energy transfer (EET) pathways: 1^-2^3 and 6 ^■ 
(5,7) — ?> 4 -> 3. Although the nearest neighbor terms 
dominate the site to site coupling, significant hopping 
matrix elements exist between more distant sites. This 
indicates that coherent transport itself may not explain 
why the excitation energy transfer is so efhcient. 

We adopt the spin-boson model to describe the inter- 
actions between excitations and surrounding protein en- 
vironments, 



^'^b = ^^ 9i3 \i) {i\{a] + aj) 



(3) 



where gij represents the coupling constant between site 
i and the mode j of the environment. The Hamil- 
tonian Hsb describes the modulation of on-site energy 
by the environment in the liner case, where the on- 
site energy on the dimensionless environment coordi- 
nate qj ~ (a J -t- aj). Such a site-environment coupling 
can model the site-vibration interaction, leading to de- 
coherence in the FMO. This decoherence can be de- 
scribed by the master equation. To derive the master 
equation, we adopt the following assumptions. (1) The 
time evolution of the whole system (the FMO -I- envi- 
ronment) is unitary, where the full Hamiltonian is as- 
sumed to be time- independent and consists of three parts, 
namely the system Hamiltonian iJ, the bath Hamilto- 
nian Hh — ^ ■ hcdja-aj and the interaction Hamiltonian 
Hsb- The deviation of the master equation is equiv- 
alent to finding the dynamics of the FMO by tracing 
out the degrees of freedom of the environment. This is 
not always possible and we shall assume that, (2) the 
system-environment interaction is sufficiently weak, so 
that perturbation theory is applicable. Moreover, we as- 
sume that, (3) the whole system is in a product initial 
state, and (4) the environment has short memory in the 
sense that the correlation time is very short. For details, 
we refer the readers to Ref . [1^, [13 • With these assump- 
tions, we can derive a master equation to describe the 
dynamics of the FMO complex, which is given by. 




FIG. 1: The excitation transfer obtained from the master 
equation. The decoherence rates are optimized for maxi- 
mal transfer efficiency, which gives 71 — 1.2,72 = 24,73 = 
0.9333,74 = 6.5333, 75 = 53.1,76 = 1.8,77 = 27 and T = 44, 
with transfer efficiency pg ~ 0.9256. The exciton is assumed 
initially on the site 1. 



where the Liouvillian takes, 
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PjPit)Pj - IPiPit) - ip(i)F, 



£38 (p) - r 



Ps3pit)P38 " ^P3p(t) - lp{t)P3 



, (5) 



dp 
'dt 



-i[H,p] + C{p)+C38{p) 



(4) 



with P, = |j)(j|, j = 1,2, 3,. ..,7, and P38 = |3)(8| = 
Pg3. C{p) describes exciton decoherence due to the site- 
environment couplings, and £38 (p) characterizes the ex- 
citation trapping at site 3 due to interactions with the re- 
action center, labeled as site 8. The presence of this 8th 
BChl chromophore has been suggested in each subunit 
of the FMO complex by recent crystallographic datajl8|. 
Furthermore, experimental data and theoretical studies 
indicated that the 8th BChl is the closest to the base- 
plate and should be the point at which energy flows into 
the FMO complexdi,!!^. 

We shall use the population pg at time T in the reac- 
tion center given by ps{T) = Tr(|8)(8|p(T)) to quantify 
the excitation transfer efficiency. Clearly, the Liouvil- 
lian Lssip) plays an essential role in the excitation trans- 
fer. Through this term the decoherence with rates "fj 
(j = 1, 2, ..., 7) can enhance the excitation transfer as we 
show in the next section. 



III. RESULTS 

It is believed that the completely coherent dynamics is 
not most ideal for the excitation transfer, and the quan- 
tum coherence itself cannot explain the very high exci- 
tation transfer efficiency. Indeed, our numerical simula- 
tions show that r can be optimized to 87.14 with 7^ = 
to obtain the excitation transfer efficiency ps = 0.6781. 

We optimize the decoherence rates T and 7^ , {j ~ 
1, 2, ..., 7) for the maximal excitation energy transfer effi- 
ciency, and find that T = 44 and 71 = 1.2,72 = 24,73 = 
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FIG. 2: The dependence of the transfer efficiency on the ini- 
tial states. Top panel: the initial state site is a superposition 
of |1> and |6). Namely, |V'(t = 0)> = cos6l|l)+sin6'exp(j</>)|6>. 
Lower panel: the initial state is \tp{t — 0)) = cosS|l) + 
sin^exp(i(;/>)|2). The decoherence rates are optimized for max- 
imal transfer efficiency, i.e., they take the same values as in 
FiglU 0.5;0.25,0.75 



FIG. 3: The transfer efficiency versus initial states: the case 
of classical superposition. The decoherence rates are the same 
as in Fig. [T] The initial state is pll)(l| + (f -p)|6){6| for red- 
dashed line, and p|l)(I| -|- (I — p)|2)(2| for blue-dotted line. 




0.9333,74 = 6.5333, 75 = 53.1,76 = 1.8,77 = 27 yield 
a transfer efficiency ps — 0.9256. With these decoher- 
ence rates, the population on each site as a function of 
time is plotted in Fig. [1] Two observations can be made 
from Fig. [TJ (1) The population shows oscillatory be- 
haviors during the first 0.5ps and then they change in a 
smooth way, (2) the excitation on site 1 and 2 dominates 
the population, while the population at site 3 is almost 
zero. The first observation can explain why the exci- 
tation transfer lasts for longer time than the quantum 
coherence (quantum oscillation), and the second obser- 
vation suggests that the site 1 and site 2 play important 
role in the excitation energy transfer. The excitation on 
site 3 is trapped and transferred to the reaction center al- 
most immediately as to obtain a high transfer efficiency, 
as Fig. [1] shows. 

Spatial and temporal relaxation of exciton shows that 
the site 1 and 6 were populated initially with larger con- 
tribution 15|. Then it is interesting to study how the 
initial states affect the excitation transfer efficiency. In 
the following, we shall shed light on this issue. Two cases 
are considered. First, we calculate the excitation trans- 
fer efficiency with exciton initially in a superposition of 
site 1 and 6, and study the effect of transfer efficiency on 
the initial states. Second, we extend this study to initial 
states where sites 1 and 2 are initially excited. These 
calculations are performed by numerically optimized the 




FIG. 4: The dependence of the transfer efficiency on the ini- 
tial states. The decoherence rates are optimized for maximal 
transfer efficiency. The upper and lower panel are for different 
initial states with (j) = 0. Upper panel: p{t = 0) = p\ipo){'4^o\ + 
(I - p)\tpo}{tpo\ with IV'o) = cos6l|l) + sin6lexp(4)|2) and 
I'tpo') — cosff exp(i0)|2) — sin 6*11). Lower panel: |2) < — > \6). 



decoherence rates for maximal transfer efficiency at time 
T — 5ps, selected results are shown in Fig. [HEIIS] 

In Fig. [2] we present the transfer efficiency as a 
function of 9 and </>, which characterize the pure ini- 
tial states of the FMO complex through \ip{t — 0)) — 
cos 6*11) + sin 9 exp{i(f))\6) (upper panel) and \ip{t = 0)) = 
cos0|l) + sin6'exp(i0)|2) (lower panel). We find that a 
properly coherent superposition of sites 1 and 6 (or 2) 
can enhance the exciton transfer. For exciton initially 
in a superposition of 1 and 6, the transfer efficiency is 
more sensitive to the population ratio (characterized by 
9) but not to the relative phase cj). The transfer effi- 



ciency arrives at its maximum with 9 



and 







(it is not unique), suggesting that exciton on site 6 is 
more favorable for the energy transfer. This finding is 
changed when the relative phase (j) is not zero, for exam- 
ple, when (/) = TT, the maximal transfer efficiency moves 
toward larger 6. For exciton initially excited on site 1 
and 2, both the population ratio and the relative phase 
affect the energy transfer, the transfer efficiency reaches 
its maximum with 6 — ^ and 6 — Q and its minimum 



with 6 ^^ and 



= 0. 



Fig. [3] shows the dependence of the transfer efficiency 
on the mixing rate p, where the initial state is a classical 
mixing of site 1 and 6 (or 2). Obviously, the population 
mixing favors the transfer efficiency. In particular, the 
classical mixing of site 1 and site 6 (or site 2) always 
increase the transfer efficiency, this is different from the 
case where the mixing is performed in a different way, 
see Fig. SI We observe from Fig. 2] that the transfer 
efficiency attains a minimal value ai 9 = jTt and p = 1 
(upper panel) which is consistent with the observation 
emerged from Fig. [51 As p decreases, the transfer effi- 
ciency increases linearly. Similar features can be found 
from the lower panel of Fig. 2] These results suggest 
that the population on site 6 favors the energy transfer: 
if the mixing of site 6 and site 1 is classical, the transfer 
efficiency increases linearly with the population on site 
6. Otherwise, it depends on the mixing angle 9 and the 
relative phase (p. 

We would like to note that the coherent superposition 
of exciton on the sites 1 and 6 (or 2) can help the energy 
transfer even when the decoherence is absent. For exam- 
ple, our numerical optimizations show that the transfer 
efficiency can be increased to 98% starting with an initial 
state cos(1.26327r)|l) -f sin(1.26327r)|2), however, this is 
not the case for classical mixing (or incoherent superpo- 
sition) of the sites 1 and 2 (or site 6) as initial states. 
The transfer efficiency is always smaller than (or equal 
to) 67.81% (the efficiency in the case with exciton 100% 
on the site 1 in the absence of decoherence) with initial 
states in the form p|l}(l| + (1 — p)|2)(2|. This suggests 
that the transfer enhancement by coherent superposition 
of exciton on site 1 and site 6 (or site 2) is a common 
feature for the energy transfer in the FMO complex re- 
gardless of the presence of decoherence. But the enhance- 
ment of transfer efficiency by classical mixing works only 
in the presence of decoherence. This can be understood 
as the coherent cancelation or coherent enhancement in 
the dynamics of exciton transport. 

It has been shown that the local decoherence can en- 
hance the transfer efficiency [ll| in the FMO, this can be 
understood as the fluctuation-induced-broadening of en- 
ergy levels, which bridges the on-site energy gap and the 
coupling between them. This gives rise to the following 
question: how the fluctuations in the site energies and 
couplings affect the transfer efficiency. In the following, 
we will study this issue and show that the exciton trans- 
fer in the FMO complex remains essentially unaffected 
in the presence of random variations in site energies and 
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FIG. 5: Effect of fluctuations in the Hamiltonian on the trans- 
fer efficiency. The upper panel is for the ffuctuations with zero 
mean, while the lower panel is for ffuctuations with positive 
mean. The resulting fidelity is an averaged result over 20 
independent runs. 




FIG. 6: Effect of fluctuations in the Hamiltonian on the trans- 
fer efficiency. The ffuctuations are Gaussian with mean /i and 
the standard deviation a. (a) The fluctuations happen only 
in the on-site energies, (b) the fluctuations in both the on-site 
energies and couplings. The resulting ffdelity is an averaged 
result over 100 independent runs. 



inter-site couplings. This strongly suggests that the ex- 
perimental results recorded for samples at low tempera- 
ture would also be observable at higher temperatures. 

We take the decoherence rates that maximize the 
transfer efficiency for the numerical simulation. Two 
types of fluctuations in the site energies and inter- 
site couplings are considered. In the first one, it has 
zero mean, while the another type of fluctuations has 
nonzero positive mean. For the fluctuations with zero 
mean, the Hamiltonian in Eq. ^ takes the follow- 
ing changes, Hjj -^ Hjj{l + ri ■ (rand(l) — 0.5)) and 



Hi^j — >■ Hi^j{l + r2 ■ (rand(l) - 0.5)). For the fluctua- 
tions with nonzero positive mean, the Hamihonian takes 
the following changes, Hjj -^ Hjj{l + ri ■ rand(l)) and 
Hi^j — !> Hi^j{l + r2 • rand(l)), where rand(l) denotes 
a random number between and 1. So a 100% static 
disorder may appear in the on-site energies and inter-site 
couplings. 

With these arrangements, we numerically calculated 
the transfer efliciency and present the results in Fig. 
m Each transfer efficiency is a result averaged over 20 
fluctuations. Two observations are obvious. (1) As the 
fluctuations with zero-mean in the site energies and cou- 
plings increase, the transfer eflficiency fluctuates greatly. 
(2) For fluctuations with nonzero positive mean, the ef- 
flciency increases with r2 but decrease with ri . This can 
be interpreted as follows. The energy gap between neigh- 
boring sites blocks the energy transfer, whereas the inter- 
site couplings that represent the overlap of different sites 
favor the transport. As a result of competition, the effi- 
ciency increases with ri but decreases with r2- For the 
fluctuation with zero mean, the results are subtle, it de- 
pends seriously on each run, and the averaged transfer 
efficiency remains almost unchanged. 

Ref.[15i] suggested that the fluctuations in the on- 
site energies are Gaussian. To examine the effect of 
Gaussian fluctuations, we introduce a Gaussian function, 

1 ("-f)^ 

y[x\a-,fi) — — y==e ~^^^~ , where ^i is the mean, while 

a denotes the standard deviation. Furthermore we as- 
sume that the fluctuations enter into only the on-site 
energies, namely, Hjj is replaced by Hjj{l -t- y(a;|cr, /x)), 
j = 1, 2, ..., 7. With these notations, we plot the effect of 
fluctuations on the transfer efficiency in Fig. HJa). We 
flnd that the efficiency is almost independent of the vari- 
ation, but it decreases as the mean increases. This again 
can be understood as the blockage of energy transfer by 
the energy gap between the neighboring sites. When 
these Gaussian fluctuations occur in both the on-site en- 
ergies and the couplings, we flnd from FiglHl^b) that the 
transfer efficiency increases as the mean increases and the 
variation decreases. This is a result of competition be- 
tween the energy gaps and couplings of the neighboring 



sites. 

Some one may wonder about the decoherence rates 
that yield the high energy transfer efficiency - why the 
optimal decoherence rate is different for each sites? why 
some decoherence rates are larger than the energy gaps 
between the sites. As a phenomenal model presented 
here, we can not answer these questions, this stimulate 
further works on this topic. 



IV. CONCLUSION 

In this paper, we study the excitation transfer effi- 
ciency in light-harvesting complexes, in particular, the 
Fenna-Matthew-Olson (FMO) complex. This is a prob- 
lem of long interest, and timely due to the recent activ- 
ity sparked by suggestions of quantum effects in energy 
transport. It is well known that photosynthetic organ- 
isms operate at extremely low light levels, and the ex- 
citon initially excited in the FMO complex may occupy 
the site 1 and 6 (or site 2) simultaneously. This stimu- 
late us to study the effects of initial state on the transfer 
efficiency. Based on the Frenkel exciton Hamiltonian and 
the master equation, we have optimized the decoherence 
rates for maximal energy transfer efficiency, which can 
reach about 93% at time 5ps. By considering different 
mixing of exciton on site 1 and site 6 (site 2) as the ini- 
tial states, we have examined the effect of initial states 
on the energy transfer efficiency. The results suggest that 
a classically mixing of site 1 and site 6 (or site 2) always 
enhance the energy transfer efficiency, whereas a coher- 
ent superposition of site 1 and site 6 (or site 2) increases 
the efficiency only for special initial states. This can be 
understood as coherent cancelation or coherent enhance- 
ment in quantum physics. 
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